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LAUNCHING OF SURFACE WAVES ON AXIAL- CYLINDRICAL REACTIVE SURFACE 


by Norman C. Wenger 
Lewis Research Center 


SUMMARY 

The excitation of the dominant transverse magnetic (TM) surface wave on an 
axial-cylindrical reactive surface is discussed. The surface wave launcher 
consists of a perfectly conducting, infinitely thin cylindrical surface coaxial 
with an axial-cylindrical reactive surface. The reactive surface is of infinite 
extent, and the perfectly conducting surface is of semi-infinite extent. The 
surface wave is excited by the dominant TM mode in the coaxial portion of the 
structure . 

Numerical results are obtained for the energy transported by the reflected 
field, the surface wave field, and the radiation field by using an exact anal- 
ysis. The exact results are then compared with the results from two approxi- 
mation techniques that are frequently used to solve problems of this class to 
determine the validity of these techniques. 

This method of excitation was found to be very efficient over a large 
range of frequencies and over wide variations in the surface reactance. 


INTRODUCTION 

The propagation of electromagnetic surface waves on various types of 
structures has been widely treated in the literature. This extensive treatment 
has been motivated largely by the inherent low attenuation and large bandwidth 
found in many surface wave structures. Excellent surveys of the properties of 
surface waves and surface wave structures useful in communications systems have 
been presented by Zucker (ref. l) and Barlow (ref. 2). A considerable amount 
of work on new structures that can support surface waves such as anisotropic 
ferrites and plasma columns is also being performed. A knowledge of the proper- 
ties of the surface waves associated with these structures can provide some in- 
sight into their composition. 

A common problem to all of these areas of interest is the efficient excit- 
ation or launching of the guided electromagnetic surface waves. A general 
requirement for a good surface wave launcher is a high launching efficiency 
over a large frequency bandwidth. Since the surface wave fields are of infinite 
extent, the launcher must also be of infinite extent to have 100 percent launch- 
ing efficiency. Brown (ref. 3) has shown that the launching efficiency of a 



finite-sized launcher can be made arbitrarily close to 100 percent. This large 
efficiency can be realized, however, only at the expense of frequency bandwidth. 


A very limited amount of work has been done in analyzing finite-sized 
launchers that are physically realizable. The majority of numerical data pres- 
ently available is for the class of launchers that are infinitesimal in some 
dimension. To this class belong the short electric and magnetic current ele- 
ments, line sources, current loops, etc. The finite-sized launcher can be 
handled, at least in theory, by a superposition of infinitesimal sources. In 
practice, it is usually difficult to carry out this superposition because of 
the complexity involved in the calculation and the uncertainty in the distribu- 
tion of the field within the launcher. The finite-sized launcher is usually 
analyzed by using an approximate field distribution. The field or aperture 
distribution is often approximated by a "chopped" surface wave distribution; 
that is, the field in the aperture plane of the launcher is assumed to have the 
same form as the surface wave field within the aperture and is assumed to 
vanish everywhere outside of the aperture. Another often-used approximation 
technique is Kirchhoff*s approximation. In this method the aperture field is 
assumed to be of the same form as the unperturbed incident field. For either 
case, the surface wave amplitude can be easily computed by an integration over 
the aperture plane since the surface wave modes and the radiation field are 
orthogonal (ref. 4). The accuracy of the results obtained by using these ap- 
proximation techniques is usually unknown since no criterion exists which can 
determine the extent of the approximat ions . 

The purpose of this report is to present an exact analysis and numerical 
results for the launching characteristics, radiation pattern, and frequency 
bandwidth of a finite-sized launcher that is physically realizable. These re- 
sults will then be compared with the results for the "chopped" surface wave 
distribution and with the results using Kirchhoff 's approximation to determine 
under what conditions the approximation techniques are valid. 

The analysis will be restricted to the class of surface wave structures 
that have the configuration of a cylindrical column of circular cross section. 
In order to keep the results of this work as general as possible, the effect of 
the surface wave structure on the electromagnetic field will be taken into 
account by specifying a surface impedance. The surface impedance is defined as 
the ratio of the tangential electric field to the component of the tangential 
magnetic field perpendicular to the electric field along the surface of the 
wave guiding structure. The numerical value of the surface impedance will, of 

course, depend on the composition 
of the structure and the polariza- 
tion of the electromagnetic field. 
The case where the structure has no 
losses will be considered to sim- 
plify the analysis; consequently, 
the surface impedance will be a 
pure reactance. Only the case of 
inductive surface reactance will be 
treated. 


The particular configuration 
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that will be considered is shown in figure 1. The structure consists of a re- 
active cylindrical surface of radius a and of infinite extent in the 
z-direction. Coaxial with this cylinder is an infinitely thin, perfectly con- 
ducting surface of radius b for z < 0. The surface wave field, radiation 
field, and the reflected field will be computed for the case where the incident 
field is the dominant transverse magnetic (TM) mode in the region a < r < b 
and z < 0, propagating in the positive z-direction. The analysis will be re- 
stricted to the frequency range where the dominant TM mode is the only propa- 
gating mode. 

Before a formal solution of the problem is carried out it is instructive 
to examine the various types of waves that exist in the different portions of 
the structure. 
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SYMBOLS 

radius of cylindrical reactive surface (see fig. l) 

radius of perfectly conducting surface (see fig. l) 

b-a (see fig. l) 

function defined in eq. (18c) 

function defined in eq. (2l) 

Weiner-Hopf factors of F(j3) 
eigenvalue for surface wave mode 
imaginary part of 
function defined in eq. (18a) 
function defined in eq. (18b) 

free -space wave number 
real part of k Q 
negative imaginary part of k Q 
eigenvalue of coaxial mode 
real part of 

radial coordinate (see fig. l) 
surface reactance 
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Z Q characteristic impedance of free space 

z axial coordinate (see fig. l) 

“ k 0 X s /Z Q 

P complex variable 

P q propagation constant for surface wave mode 

P^ real part of p Q 

p^ negative imaginary part of P Q 

X propagation constant for coaxial mode 

r' real part of y 

o o 

negative imaginary part of Y" 0 
6 azimuthal angle (see fig. l) 

p radial coordinate (see fig. l) 

cp compliment of polar angle (see fig. l) 

cp(r.,p) Fourier transform of ^(r^z) 

cp ± (r,p) single-sided Fourier transforms of i|r(r,z) 

iKr,z) 0 component of total magnetic field 

\|/^(r,z) incident field 

\|r s (r,z) scattered field 

cd angular frequency 

Subscripts : 
i inc ident 

rad radiated 

rf reflected 

sw surface wave 
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SURFACE WAVE LAUNCHER 


Coaxial Portion of Structure 

The total field in the coaxial portion of the structure can be expressed 
in the form of an infinite summation of modes. These modes will consists, in 
general, of transverse magnetic (TM), transverse electric (TE), and hybrid 
modes. The characteristics of these modes will depend on the variation of the 
surface reactance with frequency and on the polarization of the field. 

At low frequencies only the TM QO mode will propagate. The subscripts 
denoting the TM^ modes are selected such that the first subscript m corre- 
sponds to the number of cyclic variations of the mode intensity with 6 and 
the second subscript n corresponds to the number of nodes of' the mode inten- 
sity with radius. The TM QO mode is composed of an electric field with axial and 
radial components and a magnetic field with an azimuthal component. All other 
field components are zero. The remaining TM modes and all the TE and hybrid 
modes are cut off at low frequencies. Thus, the total field in the coaxial 
portion of the structure will consist of the TM QO mode except in the vicinity 
of the discontinuity at z = 0. The presence of the discontinuity will cause 
the dominant TMq 0 mode to excite the evanescent TMon modes. All the other 
possible modes are not excited since both the TM 00 mode and the discontinuity 
are circularly symmetric. The total field in the coaxial portion of the struc- 
ture will consist, therefore, of the entire spectrum of TM Qn modes with each 
mode having an axial and radial component of electric field and an azimuthal 
component of magnetic field. Since the electric and magnetic fields are re- 
lated by Maxwell ! s equations, only the magnetic field need be determined to 
uniquely specify the total field. Thus, the azimuthal or 6 component of the 
magnetic field can be regarded as a scalar function from which all the other 
field components can be derived. 

Let ^(r,z) be the 0 component of the magnetic field. Then, in the co- 
axial portion of the structure, \[r( r,z) can be expanded in a series of the 
eigenfunctions (ref. 5, pp. 709 to 778): 

>|r(r,z) = A 0 [j 1 (-jp 0 r)H^ ) (-jp 0 b) - J Q ( -jp 0 b ) h|( - jp Q r )] e 


+ 


E *n[ J l(p n r ) H o(Pn b ) " J o( p n b ) H i(£n r )] e ^ 


n=l 


( 1 ) 


where the A n are complex amplitude constants. A time dependence of e 
has been assumed for all field quantities. Since the solution (l) must satisfy 
Maxwell’s equations, the eigenvalues p n and the propagation constants y n 
are related by 


• o ^o 


+ K 


rj = p 2 


n 


k o 


for n > 0 


5 


•where k Q is the free-space wave number. The eigenvalues p n are determined 
by the boundary conditions. The ratio of the axial component of the electric 
field to the azimuthal component of the magnetic field must equal the surface 
reactance jX s along the reactive surface r = a. This condition requires 
\|r(r,z) to satisfy the equation 

JT ^ [ri|r(r,z)] + <u|r(r,z) | =° (2) 

r~a 


where a = k 0 X s /Z 0 and Z Q is the characteristic impedance of free space. 
Substituting the general solution (l) into equation (2) generates the following 
set of equations that must be satisfied by the eigenvalues p n : 


JP G a 


J 0 ( -jPo b)H g ( -jPo a) - J 0 (-jP o a)Hg(-jp o b) 

J l("jPo a ) H o("J p o b ) “ J 0 (-J p o b ) H l(-J p o a ) 


-aa 


(3a) 


N 0 (p n b)J 0 (p n a) - J 0 (p n b)N 0 (p n a) 

^ J o t Pn b N 0 Cp n b ) Ji(Pn a T 


for n > 0 


(3b) 


Along the surface of the perfect conductor r = b, the axial component of the 
electric field must vanish. This condition requires i|r(r,z) to satisfy the 
equation 

JT ^ [rt(r,z)] = 0 (4) 

r=b 

This boundary condition has been built into solution (l). 

If the surface reactance vanishes, the first term in solution (l) will re- 
duce to the ordinary TEM wave associated with a coaxial line; whereas the re- 
maining terms in solution (l) will reduce to the TM^ modes for a coaxial line. 
It should also be noted that since the total field consists entirely of TM 
modes it is only necessary to specify the surface reactance for one polarization 
of the electric field. The value of the surface reactance for other polariza- 
tions is arbitrary. 


Open Portion of Structure 

The total field in the open portion of the structure consists of discrete 
modes guided by the reactive surface plus the radiation field. Tne guided 
waves can be described as either TM, TE, or hybrid. Since the field in the 
coaxial portion of the structure is composed entirely of circularly symmetric 
TM modes and since the structure is circularly symmetric, the field in the open 
portion of the structure will also consist of circularly symmetric TM modes. 

The only transverse magnetic mode that possesses circular symmetry is the TMq 
mode, which is commonly called the Goubau wave. The single subscript n de- 
noting the various TM n modes refers to the number of cyclic variations of the 
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mode intensity "with 0. A second subscript is not necessary since all modes in 
the open portion of the structure are evanescent in the radial direction. 

The total field ^(r,z) in the open portion of the structure consists, 
therefore, of the radiation field plus the TMq surface wave (ref. 2, pp. 60 
to 69) of the form 

B o H2(-jh o r)e‘ jP ° Z (5) 

where B 0 is a complex amplitude constant. The eigenvalue h Q and the propa- 
gation constant p Q are related by 

p2 = h 2 + k 2 
o o o 

since equation (5) must satisfy Maxwell's equations. The eigenvalue h Q is 
determined by the boundary condition (2) at the reactive surface. Substituting 
equation (5) into equation (2) requires h Q to satisfy the equation 


0h o a ” 


H^( -jh a) 


H|(-jh o a) 


FORMAL SOLUTION OF PROBLEM 
Statement of Problem 

The formal solution of the problem will be carried out by using Laplace 
transform and Wiener-Hopf (refs. 6 and 7) techniques. It is convenient to 
decompose' the total field r,z) into two parts: an incident field t-^r^z) 

and a scattered field ty s (r,z), where 

t(r,z) = q(r,z) + ^(r^z) 

The incident field is the dominant TMq 0 mode in the coaxial portion of the 
structure and exists by definition for a < r < b and all values of z : 

t ± (r,z) = [j 1 (-jp 0 r)Hg(-jp 0 b) - J 0 ( -jp 0 b ) h|( - jp Q r j] e JT ° Z (7) 

The amplitude factor in equation (7) has been selected such that the incident 
field has unit amplitude at r = b. Since the incident field does not satisfy 
the proper boundary conditions for z > 0, the scattered field will contain a 
term of the same form as the incident field for z > 0 to nullify this improper 
solution. 

The scattered field satisfies the following equations : 
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( 8 ) 


a 2 t s 2. dt s d 2 t 

r c 

1 d 


3 r 2 


+ — -r — — + — —— + (k 2 _ -i- ^ 

r or ^ z 2 \ o r 2 / 


t = 0 


- ^ (rt 8 ) + 
1 d 


= 0 


r=a 


r SF (r ^s) 


= 0 


r=b 

z<0 


- 


= e 


-jr 0 2 


r=b + 

z>0 


r=b‘ 

z>0 


— ( r \|r ) 

r dr ^ 


?5W.) 


r=b + 

zX) 


= 0 


r=b“ 

z>0 


(9) 

( 10 ) 


( 11 ) 


( 12 ) 


Equation (8) is the Helmholtz equation expressed in cylindrical coordinates. 
Equations (9) and (10) are a statement of the boundary conditions on the reac- 
tive surface and on the perfectly conducting surface, respectively. Equation 
(ll) requires the scattered field to be discontinuous at r = b and z > 0 
in order to make the total field continuous. Equation (12) requires the axial 
component, of the electric field to be continuous at r = b and z > 0. In 
addition to the above boundary conditions, a radiation condition must also be 
imposed. The radiation condition requires all admissible solutions for the 
scattered field to correspond to divergent waves at infinity. 


Laplace Transformation 
Let the function cp(r,p) be defined by 

cp(r,|3) = cp + (r,p) + cp“(r,(B) 


where 


9 + (r,p) 



'|' s (r,z)e“P z 


dz 
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and 


<P"(r,p) 



^ S (r,z)e - P z dz 


The scattered field ^ s (r,z) can then be recovered by the inversion integral 


^ s (r,z) = 


sfe f 

J n 


p)eP z dp 


where C denotes a suitable contour in the complex P-plane, 


In order to make cp + (r,(3) and cp”(r,0) analytic functions of P in a 
common region in the complex p -plane , the free-space wave number k Q will be 
made complex. This is equivalent to introducing losses into the medium sur- 
rounding the structure. Let jk Q = jk^ + k Q , where k^ and k Q are real. 

In the final solution k^ will be set equal to zero to recover the result for 
the lossless case. Since k Q is complex, the propagation constants T 0 and - 
P Q will also be complex. Thus, y 0 and P Q will, in general, be given by 

jr 0 = 3*0 + *o and JPo = JPi + where *o> *o> Pi' and Pi are a11 real * 
The inequalities k^ < P^ < Tq and *o < Po - k o can be shown to be valid 
for all cases . 

Taking the Laplace transform of equations (8), (9), and (10) gives 


d2<£ 

5r^ 

+ r 5? + (p 2 + k o - = 0 

(13) 


^ (rep) + oep | = 0 

(14) 


r=a 



^lr(rcp-) 1 =0 

(15) 


r=b 


where p by virtue of the previous discussion, must be restricted to the range 

-To < Re p < r”. 


Solution for Transformed Scattered Field 
The solution of equation (13) is of the form 

cp(r,p) = A(p)J 1 (7or) + B(P)hJ(Xt) 
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1/2 

where A=(k^ + p2) / and A(p) and B(0) are suitable functions of 0 that 
must be selected to satisfy the boundary conditions (14) and (15) and the radi- 
ation condition. 


In the region a < r < b the proper solution of equation (13) is 


cp(r,p) = cp(b~,p) 


[AH^(Aa) + allf(Aa) 

? 

1 • 

1 

[AJ 0 ( Aa ) + aJj_(Aa) 

H^(Ar) 

[aHq(As) + aHj(Aa) 

? 

& i 

T i 

! 

^AJq(As) + aJ^CAa) 

H^(Ab) 


(16) 


The solution given by equation (16) is an even function of A; consequently, 
either branch of A can be selected. In the region r > b the proper solu- 
tion of equation (13) is 

H?( Ar ) 

cp(r,p) = cp(b + ,p ) -f (17) 

H|( Ab ) 

In order to satisfy the radiation condition, the branch of A where Im A < 0 
must be selected. 

The unknown coefficients cp(b~,0) and cp(b + ,p) in equations (16) and (17) 
can be determined by the discontinuity condition on the scattered field at 
r = b as given by equation (ll). This can easily be accomplished by intro- 
ducing the functions J+(b,p), J-(b,P) and E+(b,p), where 


CO 


J + (b,p) s / 

J Q 

0 s (b + ,z) 

- t s (b",z)] 

|e"P z dz 

(18a) 

J-(b,P) = f | 

-oo 

[i s (b + ,z) 

- ^ s (b~,z )] 

e”P z dz 

(18b) 

E + (b,S) = 

(rep) 


(18c) 


r=b 


It should be noted that the function J“(b,p) is actually the Laplace transform 
of the portion of the electric current on the perfectly conducting surface 
r = b associated with the scattered field. The function E + (b,p) is propor- 
tional to the transform of the axial component of the total electric field 
evaluated at r = b. Thus, E + (b,0) could have been defined equally well by the 
equation 


E+ (b,P) = ^ J— (rcp + ) 

r=b 
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by virtue of equation (15). The function J + (b,p) can be calculated at once for 
for the region Re p > -y" by using equation (ll): 


J + (b,p) 


1 

3 + jr Q 


(19) 


The functions J"(b,p) and E + (b,p) are, at present, unknown; however, a con- 
sideration of the complex forms of k Q , p Q , and y Q reveals that J"(b,p) is 
an analytic function of the complex variable p in the region Re p < f)) and 
that E + (b,p) is analytic in the region Re p > -Tq* 

From the definitions of <J + (b,p) and J“(b,p) it is apparent that 

9(b + ,P) - <p(b-,p)= J + (b,P ) + J-(b,P) (20) 

The unknown coefficients cp(b“,p) and cp(b + ,P) can be expressed in terms of the 
function E + (b,p) by using the solutions (16) and (17) in conjunction with the 
definition for E + (b,p) as given by equation (18c). The result of this opera- 
tion allows the left side of equation (20) to be put into the form 


cp(b + ,p) - cp(b“,P ) = F(p)E + (b,p) 


where F( P ) is given by 

t2/ 


*„)-j a_5<2sL 

*bA 2 H 2 (Ab) 


A + a 


h|( A a) 


H 2 (Aa) 


( 21 ) 


X A[H 2 (Aa)J Q (Ab) - H 2 (Ab)J Q (Aa)] + c^[h 2 ( Aa)J Q ( Ab) - H 2 ( Ab )<!-,_( Aa)] 

The function F(p) can be shown to be analytic in the strip -y^ < Re p < y£. 
Substituting equation (2l) into equation (20) and using the result for 
J + (b,p) from equation (19) give 


F(p)E + (b,p) = J“(b,p) + (22) 

Wi ener-Hopf factorization . - Equation (22) can be solved for the functions 
E + (b,pT and J"(b,3) by performing a Wiener-Hopf factorization on the function 
F(p); that is, F(p) will be expressed as the ratio of the two functions F^p) 
and F“(p), where F^(P) is analytic and nonzero for Re p > -y” and F“(p) is 
analytic and nonzero for Re P < y”. The details of this factorization are 
given in the appendix. It is sufficient, at this point to list some of the 
properties of the functions F+(P) and F“(P): 

(l) The function F + (p) has a single zero at P = -j3 0 and a branch point 
at 3 = - jk Q . 
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(2) The function F“(0) has an infinite number of zeros in the complex 
0 -plane and a branch point at 0 = jk D . The zeros of F _ (p) are located at 

0 = jr G and P = T n (n > 0). 

(3) The function F f (0) is of the order 0 -1 / 2 at infinity, and the 
function F~(0) is of the order 0 1 / 2 at infinity. 

The decomposition of the function F(0) into the ratio of F + (0) to F"(p) 
allows equation (22) to be put into the form 


F + (0 )E + (b,0 ) - — + jr 


F"(-jr n ) F“(0) - F~(-jr 0 ) 


P + Jr f 


+ F-(p)J"(b,0) 


Figure 2 shows the regions in the complex 0 -plane where the various transforms 
are analytic. A study of figure 2 reveals that the left side of equation (23) 
is analytic for Re 0 > -Tq and that the ri S ht slde is analytic for Re 0 < Tq- 
The equality in equation (23) holds only in the strip -Tq < Re 0 < Tq* 

Edge conditions . - The solution for the scattered field is not unique un- 
less the edge conditions are specified at r = b and z = 0 (ref. 7, pp. 75 
to 76 ) . These conditions require the axial component of the electric field to 
be of the order of z -d / 2 at the edge, which makes the transform of the elec- 
tric field (i.e., E + (b,p)) of the order of 0" 1 / 2 as 0 -* °°. A similar condi- 
tion exists for the asymptotic form of the current at the edge. This condition 
requires J“(b,p) to be of the order of p - "*" as 0 -» -°°. 


The asymptotic forms of F f (P), F"(P), E+(b,0), and J“(b,0) for large 
values of 0 show that each side of equation (23) approaches zero as 0 goes 
to infinity in the proper half plane. A function that is analytic everywhere 
in the complex 0-plane can be defined from equation (23). This function is 

equal to the left side of equation 
(23) for Re p > -Tq, the right side of 



equation (23) for Re 0 < y”, and 
either side of equation (23J in the 
strip -Tq < Re 0 < y”. Liouville's 
theorem (ref. 5, pp. 381 to 382) re- 
quires this function to be zero since 
zero is the only function that is 
analytic everywhere in the complex 
0-plane and vanishes at infinity. 
Setting the left side of equation 
(23) equal to zero gives 


F~(-jr ) 

E + (b,0) 2 

f (p)(p + jrJ 


Note that the proper edge conditions 
for E + (h (? p) are satisfied by equa- 


Figure 2. - Regions in complex (5-plane where transforms are analytic. 


tion (24). 
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Transformed scattered field . - The coefficients cp(b“^p) and cp(b + ,3) can 
now be expressed in terms of the known function E + (b,P) by using equations 
( 16 ) , ( 17 ) , and ( 18c ) : 


cp(b",p) = 


H |(V k o + P 2 b ) F ~(-jT Q ) 


(k^ + p 2 ) 1 ^ 2 Hg^kg + p 2 b)F*-(p)(p + jT 0 ) 


p"( - jr 0 ) 

F“ ( p ) ( p + jy 0 ) 


(25) 


. _ r (-Jr 0 ) 4(Vk| + 7 t) 

> P) - 1/2 / / \ 

(kg + p 2 ) F + (P) H o(^ k o + P 2 k )(P + Jr Q ) 


(26) 


Combining the expressions for cp(b - ,p) and cp(b + ,3) with the solutions (16) 
and (17), respectively, completes the solution for the transformed scattered 
magnetic field. 


Scattered Field 


Inversion integral . - The scattered magnetic field ty s (r,z) can now be 
computed by using the inversion integral 


^s( r J z ) 



cp(r,p)eP z dP 


The inversion contour C must be lo- 
cated in the strip < Re p < Tq ? 

as shown in figure 3, and be on the 
sheet of the Riemann surface that 
corresponds to the choice 
Im(k2 + p2)l/2 < o in this strip. 

The branch cuts were selected as 
straight line segments extending ra- 
dially from the branch points 
0 - ±jk 0 . For the lossless case the 

D ~ branch cuts will be located on the 
p . n 

imaginary p -axis . 

Discrete spectrum . - In the com- 
putation of the discrete portion of 
the spectrum, the amplitudes of the 
modes need be evaluated at only one 
radius since the radial variation of 
these modes is already known. It is 
convenient to evaluate the field at 
Figure 3. - Inversion contour in complex p-plane. the radius r = b since the radially 
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varying factor in the general solutions (16) and (17) has been normalized to 
unity at r = b. 

First, consider the region a < r < b and z < 0. The amplitude of the 
scattered field at r = b is 


J C 

where cp(b~,3) is given by equation (25). The contour C can be closed in the 
right half 3 -plane with a semicircle of infinite radius that is deformed 
around the branch cut as shown in figure 3. The integration along the semi- 
circle and along the branch cut can be shown to be zero for z < 0. Thus, the 
integral over the original contour C must equal -2 jcj times the sum of the 
residues of the enclosed poles in the right half g-plane. All of the poles of 
cp(b“,3) in the right half 3-plane are due to the zeros of F~(p). Thus, for 
z < 0 


t„(b,z) = 


F“(-jr 0 ) 


P=jr r 



Jn z 


(27) 


The first term in equation (27) represents the dominant TMq 0 mode in the 
reflected field, and the summation represents the evanescent TMo n modes in 
the reflected field. In the region a < r < b and z « 0 the reflected 
field consists entirely of the first term, the TMq 0 mode. 


The field in the region z > 0 evaluated at r = b can be found by using 
the same inversion integral. This time, however, the contour is closed in the 
left half (3-plane with a semicircle of infinite radius that is deformed around 
the branch cut. The integration along the semicircle can be shown to be zero 
for z > 0. Thus, the integral over the original contour C must equal Znj 
times the sum of the residues of the enclosed poles in the left half (3 -plane 
plus the branch cut integral. The poles of cp(b“,p) in the left half |3-plane 
are due to the zero of F + ( p) at p = -jp G and the pole appearing explicitly 
in equation (25) at p = -jr 0 * Thus, for z > 0 


t g (b,z) 


-jH|(-jh o b)F"(-jy o ) 


(k2 - p2)l7 2 H g(-jh 0 b)(r c 


3r> ) 


3 ~~<5 3o 


-JPo z 


-Jr Q z 


+ Branch cut integral 


(28) 
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The first term in equation (28) represents the TMq surface wave mode in the 
open portion of the structure. The second term is the portion of the scattered 
field that nullifies the incident field in the region z > 0 . 


Continuous spectrum . - The branch cut integral or continuous portion of 
the spectrum gives rise to the waves associated with the radiation field. The 
pattern of the far zone radiation field will be determined by the method of 
saddle point integration. In the region r > b the scattered field is given 
by 


^ s (r,z) 


1 

2 it j 



<p(b + ^) 


H |(|/kf + 02 r) e(3a 

H 1 (v^l + p2 b ) 


dp 


where cp(b + , p) is given by equation (26). At this point, it is convenient to 
map the complex p-plane into the complex v -plane with the mapping function 


P = -jk Q sin v 

with v = a + Jt] and to transform r and z into their spherical coordinate 
equivalents with r = p cos cp and z = p sin cp (see fig. l). The complex 
v-plane is shown in figure 4, where the entire p-plane is mapped into a strip 
of width jt in the v-plane. The quantities to Q 4 refer to the images 

of the quadrants of the p-plane in the v-plane. With this transformation the 
expression for the scattered magnetic field becomes 


'LCpjT) =- 


2nj 


/ 


TT p,, / . x -jk n p sin cp sin v 

H~(k p cos cp cos v)F (-jT Q )e 0 

H^(k b cos v)F + (-jk sin v)(y - k sin v) 

0 o vu o ' x ' o o / 


dv (29) 


If the observation angle cp is selected such that cos cp f 0, then for 

?( k oF 


k Q p » 1 the Hankel function H^(k o p cos cp cos v) in equation (29) can be re- 


placed by its asymptotic form 

O 

H£(k Q p cos cp cos v) ^ 


t 


k 0 pit cos cp cos v 


^-j^p cos cp cos v - 


and lower have been neglected. The 


where terms of the order of (k 0 p )~^/ 2 
scattered field is then given by 

l/ 2 -j(k oP cos(v - 9) - x) 

3 J-/ U \ / 


-/ 


^ r o )e 

s 2lt «5 / H 2 (k. b cos v)F + (-jk^sin v )(v^ - k^ sin v )(k o pir cos cp cos v )^/ 2 


o o 


dv 

(30) 
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Figure 5. - Steepest descent contour in complex v -plane. 


for » 1. The integral will be evaluated by deforming the contour C 

into the contour of steepest descent. The exponential term in equation (30) 
contains a saddle point where sin(v - cp) = 0 or where v = cp. In the vicinity 
of the saddle point, the quantity cos(v - cp) can be expanded in a Taylor series 
to give 

C 2 

cos(v - cp) = 1 - sin £ cos £ . . . 


where £e^ = v- cp. The imaginary part of cos(v - cp) will have its greatest 
rate of change along the contour that passes through the saddle point v = cp 
at an angle of ot/4 with the a axis and, in general, satisfies the equation 
Re[k 0 cos(v - cp)] = Re k Q . This steepest descent contour (SDC) is shown in 
figure 5. The exponential term in equation (30) has unit modulus at the saddle 

-(ko/2)f 2 

point v = cp and falls off as e oP ' ^ with distance £ along the 

steepest descent contour. If k Q p » 1, the dominant portion of the integral 
for ^ g (p,cp) will be the integrand evaluated at v = cp times the integral of 
the Gaussian term with respect to p. Performing the integration gives 


+«,(p/P) = 


F-(-o T Je' J(k o (5+ ’"' 2) 


itF+t -jk 0 sin cp)Hg(k 0 b cos cp ) (r 0 - sin cp)(k Q p cos cp) 


(31) 


If more terms were retained in the asymptotic expansion of the Hankel 
function and in the expansion of the exponent in equation (28), 1 l / s (p,cp) would 
contain additional terms of the order of (k 0 p) -3 /2 and lower. These additional 
terms would give a better approximation to the field, especially if k Q p is not 
extremely large, but they would not contribute to any net radiated power. 

In deforming the contour C into the steepest descent contour, some of 
the poles of the integrand may be crossed as shown in figure 5. This usually 
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happens if the observation angle is near it/ 2 . It simply indicates that near 
the reactive surface the total field also consists of the discrete surface wave 
modes . 


This completes the formal solution for the scattered magnetic field. The 
total field can now be found by adding the scattered field to the incident 
field as given by equation ( 7 ) . 


ENERGY TRANSPORTED BY FIELD 

The energy transported by the field will be computed by using the complex 
Poynting vector theorem. The only components of the total field that transport 
energy are the TM^ incident field (eq. (7)), the TMp 0 reflected field (the 
first term in eq. (27)), the TMq surface wave field (the first term in 
eq. (28)), and the far zone radiation field (eq. (3l)). 

Since the function t(r, z) is_the 6 component of the total magnetic 
field, the electromagnetic field E,H is given by 


E = 1 V xtao 

J^ e o 


H = 


Thus, the power P 


associated with the total field is 


P = 



Ixl* • d! 



(Vxta e )xt*a e 


as 


where S denotes the surface of suitable cross section of the surface wave 
structure and 'I'* denotes the complex conjugate of \|/. 

For the case of the incident, reflected, and surface wave fields, the only 
component of the curl of "^a^ that is of interest is the radial component 
given by -d\|r/dz. Since 


- if = ^oh 


a>* rf 


-OYo^rf 
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hi f 

T sw 

hz ~ ^o^sv 


the respective powers are 


*r G f 

p i - sr / ih 

° J a 


r dr 


■rf 


f 

“ £ o J 


t rf r r dr 


■ sw 


_ ffo / 

“ £ o J 


\ ij/ “ r dr 

I r SW U - L 


where the subscripts i, rf, and sw refer to the incident, reflected, and 
surface wave portions of the total field, respectively. For the case of the 
far zone radiation field, the component of the curl of that is of interest 

id., 0 

is the cp component given by — ^ (pf). Actually, only the portion of this 
component that varies as p~l is of interest. Thus, 

p ^ ^ “ “‘W 
for k G p » 1. The radiated power is, therefore. 



^radj 2 P 2 cos <P 


NUMERICAL RESULTS 

The powers associated with the various portions of the field were computed 
on an IBM 7094. Some typical results are shown in figures 6 to 8 for b/a= 2.3 
and in figures 9 to 11 for b/a = 10.0. The values of 2.3 and 10.0 for b/a 
correspond to characteristic impedances of 50 and 138 ohms, respectively, for a 
coaxial line with an air dielectric. The ratios of surface wave power, re- 
flected power, and radiated power to incident power are shown as functions of 
k Q d where d = b - a. The range of k Q d is restricted so that only the TM 00 
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Figure 7. - Percent of incident power in reflected wave against k Q d 
for constant values of surface reactance. Radius ratio, b/a = 2.3. 


Surface 

reactance, 

x s 



Figure 8. - Percent of incident power in radiated wave against k Q d 
for constant values of surface reactance. Radius ratio, b/a * 2.3. 



Incident power in surface 

Incident power in radiation field, percent Incident power in reflected wave, percent wave, percent 


Surface 

reactance, 

x s 


1.0Z, 



Figure 9. - Percent of incident power in surface wave against k Q d 
for constant values of surface reactance. Radius ratio, b/a = 10.0. 
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Figure 10. - Percent of incident power in reflected wave against k Q d 
for constant values of surface reactance. Radius ratio, b/a = 10.0. 



Figure 11. - Percentage of incident power in radiated wave against 
k 0 d for constant values of surface reactance. Radius ratio, 
b/a = 10.0. 



Figure 12. - Radiation pattern. Surface reactance, X s = 1. 0 radius ratio, b/a = 2.3. 

mode will propagate in the coaxial portion of the structure. The results show 
that this structure is very efficient in launching surface waves , even when the 
surface reactance is quite low, if k Q d is not too small. A comparison of fig- 
ures 6 and 9 shows that the larger value of b/a yields a higher launching 
efficiency for a given value of k Q d. 

The launcher is also very broad banded as evidenced by the small reflected 
power over a large range of k G d. It should be noted that the curves shown in 
figures 6 to 11 are for constant values of surface reactance. In practice, the 
surface reactance will be a function of frequency so that the actual bandwidth 
of the launcher will not be known until the surface wave structure is specified. 

The radiation pattern of the far zone field is shown in figure 12 for the 
case where b/a =2.3 and X g = Z Q . The curves have been normalized by set- 
ting the maximum value of the power density in the forward direction equal to 1. 
The radiation patterns for other values of surface reactance are quite similar 
to those shown in figure 12 . When surface reactance was increased, the beam 
width was found to become slightly smaller for a fixed value of k Q d. 

Since the energy transported by the total field is conserved, the power 
associated with the incident wave must always equal the sum of the powers asso- 
ciated with the reflected wave, the surface wave, and the radiation field. 

This fact was used to check the numerical results . 


MUNCHING E FFIC IENC Y- APPROXIMATION TECHNIQUES 

A common method for determining the quality of a surface wave launcher is 
to compute its launching efficiency as a function of frequency. The launching 
efficiency is defined as the ratio of the surface wave power to the total power 
radiated from the aperture of the launcher. Before presenting numerical values 
for the launching efficiency by using the results from the exact analysis, it is 
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instructive to consider two approximation techniques* This problem provides an 
excellent opportunity to check the validity of these techniques since an exact 
result has already been obtained. 

The first technique that -will be considered is to approximate the aperture 
distribution of the launcher with a chopped surface wave distribution. In the 
plane z = 0 the total field ^(r,o) is assumed to be given by 

'Kr^o) = H^(-jh Q r) for a < r < b 
^(r,o) = 0 for r > b 

The total power radiated from the aperture is, therefore, 

total = £T f | H i(-J h o r )| 2r ** 

° *4 


The surface wave field that is excited by this aperture distribution is of the 
form 

B 0 Hi(-oh 0 r)e' jP ° Z 

The amplitude of the surface wave B 0 can be computed by using the orthogo- 
nal properties of the surface wave modes (ref. 4). 

b 

| h£( - jh Q r ) 1 2 r dr 

- 

|H 2 (-jh o r)| 2 r dr 


The definition of launching efficiency shows that the efficiency of this aper- 
ture distribution is equal to B 0 . Since the launching efficiency is equal to 
the ratio of two integrals with each integral having the same integrand and 
same lower limit, the launching efficiency can be made arbitrarily close to 
100 percent by increasing the upper limit b of the integral in the numerator. 
This illustrates the requirement of an infinitely large aperture before 100 per- 
cent efficiency can be obtained. 

A second approximation technique that will be considered is Kirchhoff's 
approximation. For this method, the aperture distribution is approximated by 
the unperturbed incident field; that is, the total field r,z) evaluated at 
the aperture plane z = 0 is approximated by 

^(r,o) = 't ; j_(r,o) for a < r < b 




22 




Figure 13. - Launching efficiency against k Q d. Surface 
reactance, X s = 0.5; radius ratio, b/a = 2.3. 


^(r,o) = 0 for r > b 

where +^(r,z) is given by equation (7). 

For this distribution, the total power 
radiated from the aperture is 


■ total 


n T 


toe. 


/ 


^ i (r^o)| 2 r dr 


The surface wave excited by this dis- 
tribution is of the form 


B H?( -jh r)e 
o l v u o ' 


•jP 0 z 


where the amplitude B Q is given by 


/ 


(r,o)H 2 ( -jh Q r) |r dr 


B = 
o 


J |H 2 (-jh 0 r)| 2 r dr 


The launching efficiency is therefore 



1 ,b 


Po 

I f i (r^o)H 2 (-jh o r)r dr 

2 


r. 


/ |H^(-jh Q r)| 2 r dr / | 

•'a 


o ) r dr 


Schwarz* s inequality requires the ratios of the integrals to be less than 
unity. Thus, the launching efficiency is bounded from above by Po/^o* 

Numerical values for the launching efficiency, showing both the exact and 
approximate, are presented in figure 13. The results show that the launcher is 
very efficient over a large range of frequencies even if the surface reactance 
is quite low. The exact results always give a higher launching efficiency than 
that predicted by any of the approximate techniques; however, the approximation 
techniques give quite accurate results for k Q d > 1. The values of k Q d when 
the approximations fail also give a large reflected power as shown in figure 7. 
The launcher would not be useful in this range unless some impedance matching 
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techniques were employed. For all cases, the chopped surface wave approxima- 
tion proved to be more accurate than Kirchhoff 's method in giving an estimate of 
of the exact launching efficiency. 


CONCLUSIONS 

Numerical results were obtained for the energy transported by the re- 
flected, the surface wave, and the radiation fields by using an exact analysis. 
The results show that the coaxial launcher is very efficient in exciting sur- 
face waves even when the surface reactance is quite low. Both approximation 
techniques gave quite accurate results as long as the frequency was sufficiently 
high. Large discrepancies between the exact and approximate results occurred 
only when the launching efficiency was low. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, October 25, 1965. 
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APPENDIX - WIENER -HOPF FACTORIZATION 


The Wiener -Hopf factorization of the function F(p) consists of finding 
two functions F + ((3) and F"(p) such that F(p) = F 4 " (p )/F"((B ), where F + (p) is 
analytic and nonzero in the region Re 0 > -y” and F~(p ) is analytic and non- 
zero in the region Re P < The formal procedure for performing this factor- 

ization is based on the Cauchy integral formula (ref. 5, pp. 987-989). 

Consider Cauchy's integral formula 


f(p) 



f(s) 
s - P 


ds 


where (3 is a point interior to C and f(s) is single-valued and analytic 
within and on C. Let the contour C be of the form shown in figure 14. The 
contributions to the integral along Cg and Cp cancel. Thus, 


f(p) 



f(s) 

S - p 


ds + 


2* j 




The integration over Cp produces a function that is analytic everywhere in- 
terior to Cp, whereas the integration over C 3 produces a function that is 
analytic everywhere exterior to C 3 . This procedure allows an arbitrary func- 
tion to be decomposed into the sum of two functions with each function being 
analytic in different but overlapping portions of the complex plane. 

The same idea can be extended to the case where f(p) is analytic in an 
infinitely long strip. If f((3)/p vanishes as |(3j -» f(p) can be decom- 
posed into the sum of two functions f + (P) and f-(p) by performing the indi- 
cated integrations over and in figure 15. 

f(p) = f + (p) + f(p) 


f~(p) 


1 

2 rtj 


f + (p) 


1 
2* j 



ds 


ds 


The function f“(p) is analytic for Re P < sp, and f + (P) is analytic for 
Re p > sg. If f(P) is an even function of "P, the functions f + (0) and f-(p) 
are related by f + (-P) = f~(p) and f + (p) = f~(-P). This relation provides a 
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Re s 



very simple method for computing f + (0). 


s 2 




C 


1 


S 1 



-♦Strip where f(s) 
is analytic 


Figure 15. - Complex s-plane. 
for example, if f~(0) is known. 


In the problem under consideration, it is necessary to decompose the func- 
tion F(P ) into the ratio of two functions F+(p) and F“(p ) rather than into 
the sum of two functions, as discussed in the preceding paragraphs. The de- 
sired decomposition is easily carried out if the function In F(p) is identified 
with f(p). Since 


In F(p) = In F*(p) - In F“(p) 


and 


f(p) = f + (p) + f-(p) 


it follows that 


F"(p) = e“ f ^ 

where f + (p) and f“(p) are computed with f(p) = In F(p) by using the pre- 
viously developed formulas. Since F(p) is an even function of P, F + (p) and 
F“(p) are related by F + (p)F“(-p) = 1 and F + (-P)F“(p) = 1. 

Because of the complexity of the function F(p) in the problem under con- 
sideration, it is convenient to consider the various portions of the function 
separately. The function F(p) will be expressed, for convenience, as the 

product of the functions K(p), L(p), M(p), N(p), and (p 2 + T§) (see 
eq.. (21)), where 


K(P) 


2j 

irtbX 
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L(P) 


Hq( Aa ) 
H 2 (Ab) 


M(p) = A + a 


H 2 (Aa) 

Hg(Aa) 


N(P) = - 
A 


a2 + v-2 

o 


2 [Ho(Aa)J 0 (Ab) - Ho(Ab)J 0 (Aa)] + Aco[i 2 (Aa)J 0 ( Ab) - H 2 (Ab) J^Aa)] 


and 

A = (k§ + P 2 ) 1 ^ 2 

Firsts consider the function L(fi) ~ L + (p)/ir(0). 
given by 


l-O) - 


The factor L"(p) is 


where 


r(p) 




Figure 16. - Contour Cj in complex s-plane. 


and 


ft = (k 2 + s 2 ) 1 / 2 

It is necessary to include a convergence 
factor in the expression for Z“(p) 

since the function L(0 ) does not have the 
proper behavior at infinity. The contour C-j 
must be located in the strip -Tq < He s < 
as shown in figure 16 with the point s = 0 
located to the left of C-^. 

The integral over Cj_ can be evaluated 
by using Cauchy ! s residue theorem. The con- 
tour can be closed in the right half s plane 
with a semicircle of infinite radius that is 
deformed around the branch cut. It can be 
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shown, after a considerable amount of work, that the integral over the semi- 
circle of infinite radius and over the circle of infinitesimal radius that en- 
closes the branch point s = jk Q is zero. Since the integrand has no poles in 
the right half g plane ( ref . 8 ) , the integral over Cp is equal to the 
branch cut integral. The branch cut integral can be simplified by introducing 
a new real variable of integration x such that s = jk 0 x, The branch cut 
corresponds to the range in x of 1 < x < oo. After some manipulations, the 
expression for L - (g) becomes 



where P.V. denotes the principal value of the integral. The factor L + (g) is 
given by L + (g) = l/lr(-g) since L(g) is an even function of g. 


Next, consider the function M(g) = M + (g )/M"(g ) . The factorization can be 
performed on this function by taking the logarithmic derivative of M(g): 


In M*(g) - In M“(g ) = In M(g) 


1 _d_ 

M+(g) d P 


Mt(p) 


M+(g) 


w ,r(|3) = mTpT W m) 


Thus 


where 



The integral over Cj_ can be evaluated by closing the contour in the 
right half (3 plane with a semicircle of infinite radius that is deformed 
around the branch cut as shown in figure 16. It can be shown that the integral 
over this semicircle vanishes. Thus, the integral over the contour Cq is 
equal to 2jtj times the sum of the residues of the enclosed poles of m(s) 
plus the branch cut integral. The integrand m(s) has a pole at s = and 
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a second pole at the branch point s = jk 0 . After some manipulation it can be 
shown that 


M“(0) = M~(o)exp 


{/ 


0 - 


2(p - 30 ' p 


^ * / 


Branch cut integral dp 


The branch cut integral can be simplified by defining a new real variable 
of integration x such that s = jk Q x as was done for the case of L~ ( 0 ) ■ 

The expression for M~(p) then becomes 


M“(p ) = M“(o)exp 



2(p " jkj 




and 

I = V* 2 - 1 

The factor M + (0) can he easily obtained from M“(p) by using the relation 
M+(p) = 1/M - ( -P ) since M(p) is an even function of p. 

The function K(p) = K + (P )/K"(P ) can be factored by inspection. 


K~(P) = (P - j ^) 1 / 2 

Again, K + (p) is given by K + (p) = l/K“(-p) since K(p) is an even function of 
of P . 


The function N(p ) = N + (p)/N~(p) can be factored by expressing it in the 
form of an infinite product. This is possible since W(p ) is an even function 

of (k2 + p 2 )^ 2 and has singularities in the form of simple poles (ref. 5, 
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pp. 382 to 385). From the infinite product expression for N(p), the function 
W(p) can be easily identified as 


N-(3) 



e Pd/nrt 


The function N + (p) is equal to l/N“(-p) since 


N(p) is an even function of p. 


The remaining term (p^ + y2) i n ^he expression for F(j 3 ) can easily be 
factored by inspection. Thus, the factor F“(p) is simply the product of the 
expressions found for K“(p), L"(P), etc., and the factor F + (p) is the product 
of the expressions for KT^p), L + (p), etc. The function F + (p) has a single 
zero at P = -jP Q , because of the zero of M^p), plus a branch point at 
p a -jk 0 , because of the branch points of KT^p), L + (p), and M*(p). The 
function F”(p) has an infinite number of zeros in the P -plane located at 
p = jy Q and p = y n (n > 0), because of the zeros of W"(p), and a branch point 
at p = jk Q , because of the branch points of K“(p), L - (P), and M“(p). 

Thus far, the function F(p) has been decomposed into the ratio of the two 
functions F + (p) and F“(p) such that F + (p) is analytic for Re p > -Tq and 
F“(p) is analytic for Re p < Tq. The functions F + (p) and F - (p) are not 
unique. Both F+(p) and F-(p) can be multiplied by any function p(p) that is 
analytic everywhere in the finite complex p -plane to generate a new set of 
functions F + (p) and F“(p). The proper function p(P ) to select is that func- 
tion that gives the functions F f (p) and F - (p) algebraic behavior at infinity 
rather than exponential behavior. This selection is necessary to ensure that 
the field satisfies the proper edge conditions. 


A very lengthy and tedious study of the asymptotic forms of F+(p) and 
F”(p) reveals that the proper function p(p) is given by 


p(p) = exp 



1&A 

2 


, P d -i | 
+ — lnl 


' k o d \ 
2 n ) 


Pd 


rgd 

n . 


where y is Euler's constant and 
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